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ABSTRACT 

We present a general analytical procedure for computing the number density of voids 
with radius above a given value within the context of gravitational formation of the 
large scale structure of the universe out of Gaussian initial conditions. To this end we 
develop an accurate (under generally satisfied conditions) extension of the uncondi- 
tional mass function to constrained environments, which allows us both to obtain the 
number density of collapsed objects of certain mass at any distance from the center 
of the void, and to derive the number density of voids defined by collapsed objects. 
We have made detailed calculations for the spherically averaged mass density and 
halo number density profiles for particular voids. We also present a formal expression 
for the number density of voids defined by galaxies of a given type and luminosity. 
This expression contains the probability for a collapsed object of certain mass to host 
a galaxy of that type and luminosity (i.e. the conditional luminosity function) as a 
function of the environmental density. We propose a procedure to infer this function, 
which may provide useful clues as to the galaxy formation process, from the observed 
void densities. 

Key words: cosmology:theory — dark matter — galaxies:statistics — large-scale 
structure of universe — methods:analytical — methods:statistical 
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1 INTRODUCTION 

It is well known that the distribution of galaxies in the Uni- 
verse is not uniform. The galaxies are distributed in fila- 
ments and clusters, leaving large regions devoid of bright 
galaxies. These regions are known as voids. 

The first gian t void was the so-called Bootes void found 
bv lKirshner et aP il98ll) . Subsequently, thanks to large red- 
shift surveys, a large amount of these regions were found and 
analyzed (e.g. de Lapparent et al. (1986) and Vogeley et al. 
(1994) in the CfA redshift survey; El-Ad et al. (1997) in the 
IRAS survey; Miiller et al. (2000) in Las Campanas Redshift 
Survey; Plionis & Basilakos (2002) in the PSC^; Croton et 
al. (2004) and Hoyle & Vogeley (2004) in the 2dFRGS; Con- 
roy et al.(2005) in the DEEP2 redshift survey) 

Many of the theoretical works in the literature about 
voids are based on cosmological N-Body simulations. The 
first simulations of the dark matter distribution (Davis et al. 
1985) quaUtatively showed the existence of large low density 
regions, but detailed studies of these regions with better res- 
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olution are just becoming available (Van de Weygaert & Van 
Kampen 1993; Gottlober et al. 2003; Colberg et al. 2004 us- 
ing N-Body simulations and Mathis & White 2002; Benson 
et al. 2003 using semi-analytical models). Aside from the 
simulations, there are many analytical works dealing with 
underdensities in the mass field (see e.g. Sheth & van de 
Weygaert 2004, and references therein) that give good de- 
scriptions of voids. On the other hand, analytical works that 
define the voids by observable objects (i.e. in point distribu- 
tions) are not common in the literature (White 1979; Otto 
et al. 1986; Betancort-Rijo 1990). 

An important point about voids is the study of their 
contents. Despite the word, voids, of course, are not empty. 
The first detections of galaxies inside voids were spirals near 
the 'border' of previously defined voids like the Bootes (Dey 
et al. 1990; Szomoru et al. 1996a and 1996b). However, ex- 
trapolating the morphology-density relation (Dressier 1980) 
one might expect a population of dwarf galaxies well in- 
side the voids. Even though such galaxies have not been ob- 
served yet, they will provide, along with the voids statistics, 
a strong test for the galaxy formation models (Peebles 2001). 
There is some ongoing progress in the studies of galaxy pop- 
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ulations in voids, thanks mainly to the contribution of recent 
large redshift surveys (see e.g. Rojas et al. 2004; Patiri et al. 
2005b). 

There are strong discussions about what exactly con- 
stitute a void and therefore how to define them. In general, 
authors define what is a void depending on the studies they 
are carrying on. In some works, voids are defined as under- 
densities in a continuous underlying field (Van de Weygaert 
& Van Kampen 1993; Aikio & Maehoenen 1998; Friedmann 
& Piran 2001; Sheth & Van de Weygaer 2004; Colberg et al. 
2004). In other works, voids are irregular regions delimited 
by some kind of galaxies, the so-called 'wall' galaxies (e.g. 
El-Ad & Piran 1997; Hoyle & Vogeley 2002; Benson et al. 
2003; Hoyle & Vogeley 2004). Even though, these definitions 
give a very good idea, for example, about the shape of the 
galaxy distribution, they do not provide a particularly pow- 
erful tool for statistical inference. For this purpose we need a 
definition which does not smears the information contained 
in the actual object distribution. Note that if we filter this 
distribution in certain scale so as to obtain a continuous 
field and use it to define voids, or classify the objects not 
by an intrinsic criteria but by a distribution dependent one 
(e.g. 'wall galaxies'), information is smeared and the ability 
to discriminate between models by comparing observations 
with predictions is diminished. This effect is similar to what 
happens in regard to binned data: the best test using binned 
data is never better and usually worst than the best test us- 
ing the raw data. 

We define voids as maximal non-overlapping spheres 
empty of objects with mass above a given one. For example, 
we could define voids as maximal spheres empty of Milky 
Way-size halos, so that even though, the voids are empty of 
these halos, we can have dwarf halos inside the voids (see 
Fig.(l)). Otto et al. (1986) and Gottlober et al. (2003) also 
use this kind of definition. 

We will focus our work mainly on rare voids (e.g. the gi- 
ants ones) because the mean number of these voids is a very 
sensitive function of the clustering properties of the objects 
that define those regions. This implies that the available 
statistics on voids along with more general statistics like the 
counts in cell moments may be used to obtain accurate infor- 
mation about those clustering properties, providing a pow- 
erful tool to discriminate between different large scale struc- 
ture formation models. To this end the following elements 
are required: first, a handy framework to compute, for given 
clustering properties, the mean number of voids defined by 
certain kind of objects as a function of radius and, second, 
a precise characterization of the clustering properties which 
is both easy to use in the framework and physically mean- 
ingful. 

The aim of this work is to provide these elements and 
assess their efficiency. We will show how they may be used 
to infer properties of the processes whereby halos become 
galaxies of certain kind from the statistics of voids defined 
by galaxies of that kind. To this end we need to express the 
void probabilities in terms of the galaxy clustering proper- 
ties. This can be done in different ways. For example, using 
all galaxy correlation functions to characterize their cluster- 
ing properties we c ould, in princ iple, obtain the correspond- 
ing void statistics l)Whitjll979l) . However, in practice, this 
procedure is not feasible. Furthermore, even if it were, the 
information obtained about the clustering properties in this 




Figure 1. Our definition of voids: maximal non-overlapping 
spheres empty of objects with mass above M (filled circles). As it 
is shown, voids can contain objects with masses smaller than M 
(open circles). Note that we do not show for clarity the smaller 
object outside the void. 

representation do not have a direct connection with galaxy 
formation processes. In the procedure we present in this 
work, we first compute the probabilities of voids defined by 
halos with masses above a given value. This can be done an- 
alytically by combining statistics with purely gravitational 
dynamics. Then, we describe the clustering properties of the 
galaxies by their relative biasing with respect to halos of the 
same mass, which may be expressed in terms of the condi- 
tional luminosity function, and obtain an expression for the 
mean number of voids defined by galaxies. For this relative 
biasing we may either use t he ex isting semi-analytic mod- 
els (for example, IMo et alj (|20o3)) using our formalism in 
a predictive way, or use this formalism in an inductive way 
to determine that biasing from the observations. This will 
be presented in a formal way in the discussion, leaving its 
applications for a future work. 

The structure of the wo rk is as follows. In sec tion 2 we 
use an existing framework teetancort-Riidll990ri that al- 
lows us to derive the number density^ of voids of a given 
radius from the probability that a randomly placed sphere 
of that radius be empty of the objects defining the void 
(the VPF). We also show in section 2 how this probabil- 
ity may be obtained by means of an expression contain- 
ing the biasing of halos with respect to mass. In section 
3 we present an extension of the unconditional mass func - 
tion (UMF) jSheth and TormenllT999l: I Jenkins et alJl2001^ 
to constrained environments, which allow us to obtain the 
conditional mass function (CMF) which, in turn, is used to 

^ Technically, the correct term is the probability density, which 
is well defined even when it changes substantially over the local 
mean distance between voids (or relevant objects). However, to 
avoid confusion we shall use instead the term: "number density" 
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quantify the biasing of halos with respect to mass. In sec- 
tion 4 we show, by combining spherical collapse with the 
CMF, how to obtain the mean density profile both for the 
mass and for the halos within and around voids. In section 
5 we apply our formalism to different cases to obtain void 
probabilities and their mass density and halo number den- 
sity profiles and compare the results with those found in 
existing numerical simulations. Finally, in section 6 we dis- 
cuss how to use our formalism to obtain the probabilities of 
voids defined by galaxies. 



2 PROBABILITIES OF VOIDS: THE 
FRAMEWORK 

The general framework for evaluating the mean number 
densities (i.e. probability densities) of struct ures like voids 
JPolitzer and Preskil]|[l983: lOtto et alJll986l) has been ap- 
plied to the standard large scale structure models (i.e. 
Gaussian initial density fluctuation growing gravitationally) . 
The framework we use here is essentially an updated and 
extended version of the one developed by iBetancort-Riid 
il99Cfl . There, the number density, Po{r) of non-overlapping 
empty (of the objects deflning the voids) spheres with radius 
r is given by: 



J'o(r) = |^^^Po(r)(l + 0((nV)-^)) 
with 



1 dlnPo(r) 



dr 



V — —TTr 
3 



and 



Po{r)= I e 



-(fiV)(l+SN(S)) 



P{5/r) dS 



(1) 



(2) 



(3) 



where Po(r) is the probability that a randomly placed sphere 
of radius r be empty (this is the so-called Void Probability 
Function, VPF), n is the mean number density of the ob- 
jects defining the void, 5n {5) is the fractional fluctuation of 
the mean number density of these objects within a sphere of 
radius r as a function of the actual fractional mass density 
fluctuation within that sphere, 5, and P{S/r) is the probabil- 
ity distribution for the values of S within a randomly chosen 
sphere of radius r. The bias of halos with respect to mass is 
contained within 5^(5); in the non-biased case (which corre- 
sponds to the low halo mass limit) we have simply 5n equal 
to 5. 

Equation Q (without the last parenthesis) is valid for 
rare events, that is, when the mean distance between voids 
is much larger than the radius, which imply: 



nV » 1 



(4) 



In fact, when k is larger than about 3.5, equation Q 
with only the zeroth order term in the last parenthesis is 
sufficiently accurate. Extending eq.Q to smaller values of 
k (i.e. obtaining the terms of order k~^) is straightforward 
(but complex), however, eq.Q (without the last parenthe- 
sis) will be enough to our purposes because the most relevant 
constraints for galaxy formation models comes from not too 
common voids (i.e. k sufficiently large). 

The exponential in eq.lO represents the probability 
that a randomly placed sphere of radius r be empty of the 



relevant objects, when the mean fractional density fluctu- 
ation within the sphere take the value 5. Multiplying this 
quantity by the probability for a randomly placed sphere of 
radius r to have an inner mean fluctuation between S and 
S + dS, P{S/r) dS, and integrating over all possible values of 
5 we obtain Po('')- 

It must be noted that the exponential in Eq.J^ gives 
correctly the probability that the sphere be empty only when 
the clustering of the relevant objects conforms to a non- 
uniform random Poissonian process (Peebles 1980). This is 
the case to a very high accuracy both when the objects are 
mass particles or halos of a given mass, on scales (as those 
of voids) which are much larger than the size of the halos. 
For galaxies, the model might not be so good. In section 6 
we show how to modify expression ^ to be valid for ob- 
jects whose clustering properties conform to any possible 
interesting model. 

For mass particles 5n ~ S, so that eq.® is particularly 
simple. From this expression it is easy to see that as the 
density of particles (n) increases, the size of the voids goes 
to zero. 

The probability distribution, P{5/r), for a 
given power spectra is given, for any value of S in 
iBetancort-Riio Sz Lopez-Corredoira I i2002h . However, since 
the evolution of large voids is well described by the spherical 
model, we may use in eq.J^ the spherical approximation to 
P{5/r), also given in that reference with an error that for 
the voids under consideration is not very relevant, although 
in some cases, when high accuracy is required, the full 
P(S/r) may be needed. So, we shall use for P(S/r): 



1 



exp 



^(r(l + i)l/3)) 



P{S/r)dS= P{5i/r)dSt 



(1 - (1 - f )(1 - (1 + 5)-i/3))-:5 U(^(i + 5)1/3)) 



IdlnoV) ^0 54^0^731^ 
3 dlnr 



lOh-^Mpc 



(6) 



where (T^(r-) is the variance of the linear density fluctuation 
with a top-hat flltering on a scale r as a function of r; Si is 
the linear value of the density fluctuations which is related 
to S by the spherical model. For Qm = 0.3 and = 0.7 we 
have: 



6 = D{Si) = 0.993[(1 - 0.607(5, - 6.5 x 10"^ x 

X (1 - e{Si) + e{Si - i.55))<5f ))-^'«« - 1] 



where 



eix) 



1 if X > 
if X ^ 



(7) 



(8) 



Using eq.Q and eq.(|7| in eq.l^J and changing the in- 
tegration variable to 5i, we have for Po{r): 



Po{r) = / e-^^ p(^Si/r) dSi 



(9) 



where 5n{Si) is the mean fractional fluctuation within r of 
the number density of the objects defining the void as a 



© 0000 RAS, MNRAS 000, 000-000 



4 Patiri et al. 



function of the linear fractional density fluctuation within 
r. Alternatively, integrating over 5 we may write: 



-nV (1+S,^(DL{S))) 



P{5) dS 



(10) 



where 



DL{5) = 



1.68647 



1.68647 - 



1.35 



1.12431 



(1+5)2/3 (1 + 5)1/2 

0.78785 



(1+5)0.58661 

This expression for DL{5) llSheth and Tormenl liool) 

corresponds to the same cosmological parameters as eq.Q 
with Sc — 1.676 and is the inverse function of eq.Q. Note 
that we write Sn{Si{S)) rather than simply Sm{S) because it 
is 5n{Si) that we may compute directly (see next section), 
while 5n{S) is obtained through the dependence of 5i on 
S. Here we shall use the first expression (Eq.Q) for Po{r) 
while the second must be used when the exact P{S), rather 
than the spherical collapse approximation, is needed. Using 
eq.Q (or eg. llOU ) in eq.Q we obtain the mean density of 
non-overlapping empty sphere of radius r, Po{r)- The num- 
ber density of voids so that the largest sphere that they can 
accommodate have radii between r and r + dr, n{r)dr, is 
related to Po{r) by: 



Po{r) = n(> r) + 2air) + n(> (1 + -^)a2r) + 

v3 



+n(> (1 + ^/3/2)aar) + (1 + 2^/2/3)a4r) + 

+6n(^ 2.738a5r) H 

n(> r) = n{r) dr 



(12) 



where Oi depends on the mean ellipticity of the voids and 
may be taken equal to 1 without losing much accuracy. But, 
in fact, in the interesting cases [k » 1), n{2r)/n{r) << 1 
, so we may write: 

Po(r-)~n(>r) ; n{r) ~ --^Po{r) (13) 

ar 

The mean number, N{r, V), of voids with radius larger 
than r within the volume of observation V is then: 



N{r,V) = V n{> r) ^ VPo{r) 



(14) 



collapsed objects nc{m, Q, q, 5i) at a distance q from a point 
such that the mean density within a sphere of radius Q (ra- 
dius of the void) centered at this point is Si (the linear value; 
S — D{5i) is the actual one). Q and q denote the Lagrangian 
radius; for Eulerian radius (i.e. present comoving radius) we 
use respectively R and r. We then have for Sns'- 



l+Sns{r 



,s,) 



+ 



iV„(m) 



Q 



g , Nc{m,Q,q,5i) q dq 



(11) 



Nu{m) = / n„{m) dm 



Nc(m,Q,q,5i) 



nc{m,Q,q,5i) dm 



(15) 
(16) 
(17) 



where nu{m) is the unconditional number density of col- 
lapsed object with mass m. The bracketed expression is the 
mean value within the sphere of radius Q of the number 
density, Nc, of collapsed objects with masses above m when 
the mean linear fractional mass density fluctuation within 
it is 5i. 

As indicated in eg. 1151 . S„3 depends, in principle, on 
m, Q and Si. However, it may be shown on general grounds 
(and we have fully checked) that ric, Nc depend on q and 
Q almost entirely through q/Q. The reason for this lays in 
on the goodness of the approximation given in expression 
1351 1. As mentioned bellow this expression, there is a small 
residual dependence on Q, but it is completely negligible 
within the relevant range of Q values (less than a factor 2). 
Then it follows from eg. 11511 that 5ns is independent of Q, 
since, using the change of variable u = q/Q, we may write 
this equation in the form: 



1 -I- 3„s{m,Si 



1 



Nn{m) 



Nc{m, u,5i) u du 



(18) 



To obtain the fluctuation in Eulerian coordinates, 5n, 
we only need to note that as the void expands the initial 
fractional halo number density further diminishes by the 
factor {I + 5): 



l + 5jv(m,<50 = {l + 5ns{m,Si)){l + 5{5i)) 



(19) 



This is the expression that we use in eq.(|5J to obtain 
Po('"). The unconditional mass function of collapsed objects. 
nu{m) is given with high accuracy by the unconditional 
Sheth & Tormen approximation (Sheth & Tormen 1999, 
2002; hereafter ST): 



3 DERIVATION OF 5n{5l) 
3.1 Steps to follow 

To derive the mean fractional number density fluctuation 
within a sphere, 5jv, as a function of the mean linear frac- 
tional density fluctuation within it, 5i, we first obtain the 
mean fractional fluctuation of the number density of proto- 
halos within the sphere before the sphere expands in comov- 
ing coordinates, S^s, as a function of 5i. 5ns may be called 
the statistical fluctuation, since it is due to the clustering 
of the protohalos in the initial conditions before they move 
with mass. To obtain 5ns{5i) we need a framework that en- 
ables us to obtain the conditional mass function (CMF) of 



n^ST{m,5c,(T{m)) = ^ - 



2\i/2 



A 



1 + 



aS'i 



(20) 



1/2 Qh Oc 



Sc da 



■ exp 



aSc 
2a2 



where A — 0.322, p = 0.3 and a = 0.707, gt stands for the 
background density, a is the rms linear mass density fluc- 
tuation and 5c is the value of Si corresponding to collapse in 
the spherical model which for the cosmological parameters 
that we use (r2m=0.3, fl\—0.7) is 1.676. Our problem now 
is to obtain a similarly accurate approximation to the con- 
ditional mass function, nc{m,q,Q,Si), so that we can use it 
in Eg.llTHl. 
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3.2 Local constrained mass function 

3.2.1 Why do we need a new approach to the CMF ? 

Expression H18^ gives the Lagrangian constrained accumu- 
lated mass function, NcL{nn,Si), averaged within a sphere 
with mean inner linear underdensity Si and radius Q (whose 
dependence on Q has been neglected) given the local La- 
grangian accumulated mass function at a distance q from 
the center of that sphere, Nc{m,u,8i): 



NcL{m,5i) = Nu{ni){l + 5„s{m,5i)) 



(21) 



Thus, the accumulated Eulerian mass function averaged 
within the sphere, NaE, which is the ordinary mass function, 
is given by: 

Na{m,Si) = NcE{m,5i) = Nu{m){l+Sr,s{m,Si)){l+5{5i)){22) 

These expressions involve an integral within the 
sphere of the local Lagrangian constrained mass func- 
tion, nc{m,q,Q,5i). Thus, this last function is necessary 
to obtain both Nc{m,Si) and through expression (I19II . 
5^(171, 5i). There are several approaches (Sheth & Tormen 
2002; Gottlober et al. 2003; Golberg et al. 2004) giving rough 
approximations and providing reasonable fitting formula for 
NcijUySi), however none of them provides directly (with- 
out fitting) an expression sufficiently accurate to allow us 
to evaluate the void number densities, which depends very 
sensitively on 5jv(m, 5i). This is due to the fact that the lo- 
cal mass function changes substantially from the center of 
the sphere (g = 0) to its boundary [q = Q). In fact, for 
any 5i, the mean value of nc{m, q,Q, Si) /nu{ni) within the 
sphere is approximately the square root of its value at the 
center and the value of this quantity at the boundary is al- 
most the cubic root of its value at the center. So, computing 
the mass function at the center instead of the mean value, 
or assuming that the interior of the sphere may be replaced 
by a homogeneous environment characterized by the mean 
properties of the actual one do not lead to sufficiently good 
results. 

There have been several attempts to derive analytically 
the CMF. Arguably, the most motivated one is that com- 
bining the excursion set formalism (Appel & Jones 1990; 
Bond et al. 1991) with ellipsoidal dynamics (Sheth & Tor- 
men 2002). However, this procedure is not appropriate to 
our purposes because, by construction, it gives the CMF at 
the center of the sphere, which, as we stated before, is quite 
different from the mean within the sphere, which is the one 
that we need. In principle, one could repeat the ST deriva- 
tion at g = for any value of g and average over the sphere 
as indicated in Eq.(9), but this implies a rather complex 
problem that can not be solved without some approxima- 
tions. Furthermore, even if the problem could be treated ex- 
actly it would provide at most a good fitting formula where 
some parameters have to be slightly modified (with respect 
to those given directly by the formalism) to match numer- 
ical simulations, as, indeed, have already been done for the 
unconstrained case (Sheth & Tormen 2002). 

In another approach (Gottlober et al. 2003), the interior 
of the sphere is treated as a homogeneous environment and 
the unconditional mass function is rescaled to it. But, leav- 
ing aside some queries about the motivation for this proce- 
dure, in fact, it disagrees substantially with the simulations 
for large masses. 



Summarizing, neither the available CMF's nor any 
other we can envisage derived from simple considerations 
can a priori be expected to give results which are sufficiently 
accurate to our purposes. Fortunately, although we can not 
directly obtain analytically a satisfactory CMF, we may an- 
alytically extend the UMF through a procedure which is, in 
practice, exact. 



3.2.2 CMF: extending the Unconditional Mass Function 

To obtain the CMF, we simply note that, as long as the local 
evolution at a conditioned point is the same as at an uncon- 
ditioned one, the conditional local mass function of collapsed 
objects may be derived from the statistical properties of the 
local linear field in the same way as the unconditional one. 

As to the validity of this assumption, three reasons may 
be advanced: 

• Given the large difference between the scale of the con- 
straint (that of the void) and the scale corresponding to the 
masses we consider, the conditional shear distribution (of 
the field filtered on the scale of those masses) can not be 
very different from the unconditional one 

• The profile of the linear density field within the void 
is not too steep. This means that the mean negative value 
of the shear in the radial direction imposed by this profile 
is rather small (it would be strictly zero for a flat profile). 
So, the departure of the local shear distribution from the 
unconditional one is smaller than implied in general by the 
first consideration 

• It must be noted that the shear distribution plays a 
secondary role with respect to the trace (of the local ve- 
locity field tensor) in determining the mass distribution of 
collapsed objects. The difference between ST and the PS 
formalisms is due to the fact that in the former, the shear 
distribution is taken into account. The difference between 
the results of both formalisms is not that large (less than 
a factor 2). So, the small change in the shear distribution 
within a void which, according to the two previous consid- 
erations, is small, should lead to a negligible error for our 
extended mass function 

That is, if the constrained field behaves "locally" as an 
isotropic uniform random Gaussian field (with locally de- 
fined mean and power spectra), or, alternatively, if the shear 
distribution of the linear velocity field is at a constrained 
point equal to that at a randomly chosen one, we may ob- 
tain the local mass function using the Unconstrained Mass 
Function (UMF) for this local field. The relevant statistical 
property is the probability distribution, P{52/Si,q,Q), for 
the linear density fiuctuation, 52, on scale Q2 (that of the 
halos considered) at a distance g from the center of a sphere 
of radius Q (the protovoid) with mean inner linear density 
fluctuation 5i (see the conceptual diagram in flgure l|2J)- 

For a Gaussian fleld we have: 



exp 



1 _^ ■^l ' 



■;t- 

-2 '12 



Pi52/5,q,Q) = 



where 



(23) 
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gI =< si >= a^{Q); al =< Sl >= o\Q2) (24) 
1 1'°° 

<y^{x) J W^ixk) dk (25) 

ai2 = (7i2{q,Q,Q2) ^ (26) 
I Sk 1^ W{Qk) WiQ2k) sin{kq) k dk 



27r2g 



W{x) 



:{sin X — X cos x ) (27) 



where | (5^ p is the hnear power spectrum of density fluc- 
tuations. Comparing with the distribution of ^2 at a ran- 
domly chosen point which is the one implicitly involved in 
the derivation of the unconditional mass function, we find: 



P{S2 



exp 



HI) 

27r(T2 



(28) 



We see that, at least for the one point statistics, the 
field S2 at a constrained point behaves like an unconstrained 
field but with a mean value proportional to Si and a modified 
power spectra. It may be shown (Rubifio et al. 2006) that 
the joint probability distribution for the field S2 at several 
neighboring points follows very closely a Gaussian multivari- 
ate with the same mean and power spectra as the one point 
distribution. 

It is then easy to realize that the conditional mass distri- 
bution could be obtained through the following sub s titutio n 
in eg. (1201 and a renormalization fsee lRubino et all i200(t )): 



nc{m, Si,q, Q) oc nusrim, S'^, a' {m)) 



Sl 



(Ti2{m, q) 



/ 2/_N 0-12 ("T-, a) \ 1/2 

a (m) = (0-2(771) 2 j 



(29) 
(30) 

(31) 



(72, (712 depends on mass through the mass-scale relation- 
ship: 



(52(771) = 



3.51 X lOii/i-^Mp 



1/3 



ft-^Mpc 



0-2(777) = o(Q2(7n)) 
ai2{q,Q,m) = 012(9,(5,(32(771)) 



(32) 

(33) 
(34) 



after this substitution, we obtain the local Lagrangian mass 
function, nc{m, 5i, q, Q), which as we advanced, in practice, 
is only a function of q/Q. 

Note that this extending procedure is not compromised 
with any particular fit to the UMF. Actually, we could use 
for example the fit proposed by Jenkins et al. (2001). For our 
purposes, however, expressions (12) is to be preferred, be- 
cause it is very accurate in the mass range we are interested 
in. 

It must be noted that for our expression for the local 
number density of collapsed objects of mass m to be valid 
this quantity must change little within a distance of the 
order of the scale corresponding to 777. However, since the 
scale of variation of the density (for any m) is on the order of 
the void radius, it is clear that this condition holds provided 
that (52 (m)/Q « 1. 




Figure 2. In order to compute the CMF evaluated at certain 
mass at a distance q from the center of a sphere of radius Q with 
mean inner linear underdensity Si, we consider the probability 
distribution at q for the mean linear underdensity (^2) within a 
sphere with radius Q2 corresponding to the mass under consider- 
ation (see text for details). 



To carry out the computation we first obtain a fit to 
012(777, g)/o^ ((5) and find that it has the form: 



ci e 



(35) 



cri2{q,Q,m) 

where ci and C2 are certain coefficients almost independent 
of mass for (52(777) « Q, and only mildly depending on Q. 
For Q between 3.3 and 6.6/7~^Mpc, which include all voids 
we will consider, C2 goes from 0.481 to 0.554. So, we may 
use the following values for ci,C2 over the whole range: 



ci = 1.3212 



C2 = 0.525 



(36) 



inserting this in expression 1291 and using 1181 we find for 
Sns{Si,m): 



l + S„s{S,,m)=A{m)e'''^"''>'' 
A(7n,) ~ 1 ; 6(7(7) ~ 6' (777)72 



(37) 
(38) 



which provides a very good fit for 5i < 0. ^4(777) and 6(777) 
are coefficients depending only on mass (for values of Q in 
the relevant range). Note that we may need different values 
of ci, C2 for voids substantially larger than those considered 
here. So, the coefficients A, b will, in general, depend on Q 
as well as on 777. The values of b'{m) corresponding to (37) 
are given in Appendix A. 



4 MASS DENSITY AND HALO NUMBER 
DENSITY PROFILES IN VOIDS 

In computing the void number densities we have used, 
among other things, our CMF and the spherical collapse. 
Here we shall have the opportunity of checking separately 
how these assumptions works in explaining the structure of 
individual voids. 

We start with the spherically averaged mass density 
profile. A given void, characterized by its radius, R, (i.e. 
that of the largest sphere it can accommodate) and the mean 
fractional mass density fluctuation within it, So, has a unique 
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profile. However, over the ensemble of all voids characterized 
by these two parameters the density profile varies. What we 
want to obtain is the mean profile over this ensemble and 
its dispersion both parameterized by 5o and R. 

For the rare voids that we shall consider the void den- 
sity profile is equal to that for a randomly chosen sphere of 
radius R with inner fractional fluctuation So. This is so be- 
cause the fact of whether or not the sphere contains objects 
(of the type defining the void) can modify the properties of 
the profile only through the value of So, which we hold by 
construction fixed. So, the mean profile within a randomly 
chosen sphere (with fractional underdensity So) is practically 
unbiased with respect to that for a void with the same radius 
and inner underdensity. With this in mind, we may obtain 
the profile by means of the probability distribution for S{r) 
(the mean fractional enclosed density fiuctuation) at a dis- 
tance r from the center of a randomly chosen sphere with 
the condition that 5{R) = So- Transforming to the initial 
conditions with the relationship of the spherical expansion 
model DL{S) ('eg. llOll ') between the actual fiuctuation, S, 
and its linear value. Si, our problem is reduced to obtaining 
the probability distribution, in the initial field, for the value 
of Si at a Lagrangian distance r(l -I- Sy^^ (which transforms 
into present Eulerian distance r) from the center of a sphere 
with Lagrangian radius R{1 + So)^^^ and mean inner linear 
fluctuation Si — DL{So)- Since the initial conditions are as- 
sumed to be Gaussian the distribution of Si at a fixed value 
of q conditioned to DL{S{R)) — DL{So) = Si is immediately 
given by: 



exp 



P{S,/q,Si)-- 

ai = (7{Q) 
(72 = o-(g) 

o-i2(g) = 



Q = i?(l + (5)'/' 
q = r{l + Sy^'' 



T5— J 



(39) 

(40) 
(41) 



1 

2^ 



W{qk) 



X WiQk) dk 



Note that this expression gives the conditional probabil- 
ity distribution for Si at a fixed q. This is not the conditional 
probability distribution for the value of Si corresponding to 
the value of S (through Si = DL(S)), at some fixed r, which 
we represent by S(r). If it were, we could obtain immediately 
the conditional distribution for S{r) using expression (41) 
and the relationship between S and Si. The correct deriva- 
tion of this distribution can be made by a simple (but te- 
dious) argument that we give in Appendix B. We find for 
the probability distribution for S (at fixed r)conditioned to 
S{R) = So: 



P{S/r,S,) = ±P^{A) 



DL(A) 



g"i2(';) 



<5i 



V2 (cr: 



^12(9) -) 1/2 



(42) 
(43) 



<5i = DL{5o) 



g = r(l + A) 



1/3 



(44) 



with it we have for the mean profile, 5{r): 



S{r) 



S P{S/r,Si) dS ■ 



(l-Pc(A)) dA-1 (45) 



This integral must extend only to A values such that 
the integrand increases monotonically with A. Calling u{A) 
the argument of erfc in eq.(44) and ea. l45t one may check 
that the solution to equation: 



u{A) = u 



(46) 



may have more than one branch. A necessary and usually 
sufficient condition for equation 14211 to be valid (see Ap- 
pendix B) is that a branch, A'^{u), monotonically increas- 
ing with u does exist. Other branches correspond to profiles 
that have experienced a large amount of shell-crossing, so 
that expression (1421 is not valid. In order to account only 
for the relevant A^(it) branch the above condition must be 
imposed upon integral II45II . 

In an alternative procedure we may lift the mentioned 
condition on eg. 14511 (the first equation) using for P{S/r, 5i) 
the absolute value of expression 1421 and dividing it into 
the integral over S of the absolute value of expression 1421 . 
which is larger than 1 when there are additional branches. 
The difference between this procedure and the former gives 
a clue as to their accuracy. They are exact only when they 
agree; otherwise none of them is exact, the former giving 
a somewhat better result. In a similar way we may obtain 
S^{r), and the profile dispersion as{r): 



as{r) = {S^r) - S\r)) 



1/2 



(47) 



As long as the dispersion of the profile is small, which 
according to the simulations of lGottlober et all (|2003|) holds 
up to 15/i~^Mpc for the 10/i~^Mpc void, the mean actual 
profile should not differ much from the transformed mean 
linear profile, which we call the most probable profile (in 
fact, it is very nearly so). Now, from eg. 1421 and (44) we 
see that the most probable profile is essentially given by the 
center of the Gaussian (in eg 1391 i.e. the mean inner profile): 



g'i2(g) 



Si 



(48) 



where q and Q are the Lagrangian radius. Transforming 
S into Si through the spherical model relationship DL{S) 
(equation llll 'l. and using: 



o'i2(g) 

with 



Siq/Q) 



1 dLna'^iQ) 
4 dLnQ 

we may write: 

DL{S{r)) = Si S(r(l + 5(r))'''VQ) 



g(r) = r{l + S{r)) 



1/3 



(49) 



(50) 



(51) 



(52) 



This equation defines implicitly the "most probable" 
profile S{r,R,So) parameterized by R and So- This profile is 
simply the initial mean profile transformed according with 
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the spherical model. So, although presented in a somewhat 
different way, the computation is the same as those found 
in the literature (see e.g. Van de Weygaert & van Kampen 
1993). 

It follows from eg. 1421 that through the following sub- 
stitution in ea. H51^ : 

S{q/Q)5^ S{q/Q)5^ ± {a\q) - S'' {q/ Q)o\Q)f'^ (53) 

we can obtain the equations for the upper (+) and lower (-) 
68 % confidence level profiles. 

The halo number density profiles may now be obtained 
by combining the mass profiles with equation 13 711 . which 
gives the fractional fiuctuation, 5ns{5i), of the halo number 
density, previously to mass motion (i.e. due to the statistical 
clustering of the protohalos) as a function of the mean value 
of 5i within the sphere. The entire fractional fluctuation, 5n, 
is given by: 

1 + 5iv(r) = (1 + 5(r))(l + 5„.(5i(r))) (54) 

In the approximation leading to equation 15111 (i.e. 
where we simply transform the mean linear profile) the 
derivation of 5jv(r) (most probable value) is particularly sim- 
ple since, in this case, for a given r there is not only unique 
Sn and 5 but also unique 5i{r). We may then write in ea. l54ll 
the 5 value given by ea. (l5H and use DL{S{r)) (expression 
11) for Si{r), that is: 



5iv(r) = (1 + 5(r))yl(m)e'''('"'(-°^(^('''»' 



(55) 



S{r) being the solution to equation 15111 for given values of r, 
R and Sq. This gives the most probable halo number density 
profile parameterized by R and Sq. It must be noted that 
the full probability distribution for Sn at a distance r from 
the center of the void may be obtained through an argument 
similar to that leading to eg. (1421 . We find: 



P{SN/r,S{R)=So) 



1 d 



2dAi 



-erfc 



DL{Am) ~ ^<5i 



^ (^2 - Ii42l)i/2 



q = r{l + D{DL{AN))) 



1/3 



(56) 



where Si = DL{An) is the solution for Si as a function of 
A]v of the equation: 



1 + il+D{Si)){l + S„s{Si)) 



(57) 



with D(Si) given by eq.O and S„s{Si) given by eg. 13711 . The 
following relationship: 



(58) 



between the local fractional fluctuation at r, S'^{r), and the 
average enclosed fluctuation within r, Jjv(r) (also valid be- 
tween S'{r) and S{r)) may be used to obtain the proflle of 
local halo number density. 



5 RESULTS 

5.1 Void counting Statistics 

In this subsection, we apply our formalism to compute the 
mean number densities of voids for several cases, and com- 



pare them with those found by iGottlober et alJ (|20o3) by 
means of numerical simulations. In order to make a di- 
rect comparison we have applied the formalism to the cases 
treated in the mentioned work. 

They carried out high resolution TV-Body simulations 
using the Adaptive Refinement Tree code (ART) of a cube 
with 80/i~^Mpc side. The total number of particles is 1024^^ 
which leads to a maximum resolution of 4 x lO^/i~^M0 per 
particle, and a minimum halo mass of 10^/i"^Mo; the cos- 
mological parameters are i},n = 0.3 and SIa = 0.7. 

They identified voids following a criteria similar to ours; 
they considered voids as maximal empty spheres in the dis- 
tribution of dark matter halos (considered as point-like ob- 
jects). In that simulations, they found that for voids defined 
by halos with mass larger than 5 x 10^^ h~^MQ , the 20 largest 
voids have radii larger than 7.49/i~^Mpc, the 10 largest have 
radii larger than 9.2/i~^Mpc, the 5 largest, 11.0/i~^Mpc and 
the 3 largest, 11.3/i~^Mpc. 

On the other hand, when the voids were defined by ha- 
los with mass larger than 10^^h~^MQ, the 20 largest voids 
have radii larger than 6.95/i~^Mpc, the 10 largest have radii 
larger than 8.81/i~^Mpc, the 5 largest, 11.95/i~'^Mpc and 
the 3 largest, 12.63/i~^Mpc. The halo number densities (n) 
were 7.44 x 10"^(/i~^Mpc)-^ and 4.08 x 10"^(/i~^Mpc)~^ 
respectively. 

The expected number of voids with radii larger than r 
within a box of size L{= 80/i~^Mpc), N{r, L), is given by: 



(59) 



N{r,L)= I V{r)n{T)dr ~V{r)Po{r) 



V{r) = {L- 2ry 



(60) 



V{r) is the available volume for the voids (for their centers) 
that, since most voids larger than r are only slightly larger 
than r and expression (15) is a good approximation, the last 
result follows. 

Po{r) is given by expression (0 with Poir) given by 
expression 0. For (1 -I- Sn) we have: 



(I + Sn) = il + 5){l + S„s) 



(61) 



where 5„s is given by eq. I|37^ with A = 1 and b{m) = 
fe'(m)/2. For b'{m) we have used the fit given in Appendix 
A. 

In table 1 we summa rize our results and com pare them 
with the results found bv IGottlober et al] ((200^. 

We have also estimated the size of the largest void ex- 
pected in the simulation box at the 90 and 68 per cent of con- 
fidence level, V{ro)Po{ro) = 0.10 and V{ro)Po{ro) = 0.32 
respectively. These results are shown in table 2 along with 
the largest voids actually found in the simulations. 

From these results we may infer that expression Q 
gives good results for values of k over 3.5. However, for 
N > 7, regardless of the values of k, our results differ sub- 
stantially from those found in the simulations. This is due 
to the fact that the simulation box is small so that it con- 
tains only seven or so underdense structures (within which 
voids are found) rare enough (|5i|/(j ^ 3) for the spherical 
collapse to be a good approximation. To obtain good pre- 
dictions for voids such that N > 7 we must use expression 
(13) rather than eq.(15) and the full expression (Betancort- 
Rijo & Lopez-Corredoira 2002) for P{S/r) must be used in 
expression (9). However, this will rarely be necessary since 
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Table 1. Voids: our results vs. simulations. Np is the Mean num- 
ber of voids predicted by our formalism within the same vol- 
ume as the simulation s one. Nsim is Number of voids found by 
iGottlober et al] j20o3). Po is the VPF and k is the rareness of 
voids (defined in eq.(4)). The larger the value of k the rarer is the 
corresponding void. We can see here that as voids become rarer, 
the analytical predictions are in better agreement with the results 
from simulations. 



Radius 

(h-lMpc) Po Np Nsrm k 



Table 2. Largest voids. We can see that the agreement between 
our predictions and results from numerical simulations is excellent 
at 90% Confidence Level. 



Mass 
/i-^Mq 


90% CL 


^max 

68% CL 


P-max 

Simulations 


5 X 10" 


13.8 


13.0 


13.0 


1 X 10^2 


14.2 


13.55 


13.97 


2 X 10^2 


16.2 


15.38 


14.4 


5 X 10^2 


20.04 


18.6 


20.03 



mass = 5 X 10 




n=0.00744 


b= 


=0.0797 


11.3 


0.00101885 


3.4 


3 


4.880 


11.0 


001 49288 

K) .\J\J l-'-tV ^fjij 


4.8 


5 


4.621 


10.8 


0.00191639 


6.0 


7 


4.471 


9.4 


0.00984866 


24.8 


10 


3.456 


7.4 


0.072409 


119.2 


20 


2.227 


mass = 1 X 10 


l^h-lMg 


n=0.00408 


b= 


=0.1084 


12.6 


0.000427927 


1.8 


3 


5.434 


11.95 


00107413 


4.1 


5 


4.945 


11.2 


0.00191895 


7.4 


7 


5.035 


8.8 


0.0423952 


70.2 


10 


2.756 


7.0 


0.193217 


170.3 


20 


1.729 


mass = 2 X 10 


12/i-1Mq 


n=0.002162 


b= 


=0.1389 


14.05 


0.00119284 


2.2 


3 


5.432 


12.8 


0.00469144 


6.9 


5 


4.432 


11.8 


0.0125963 


15.5 


7 


3.710 


10.8 


0.0307436 


30.7 


10 


3.050 


7.5 


0.296388 


100.5 


20 


1.332 


mass = 5 X 10 


l^h-lMg 


n=0.000922 


b= 


=0.3010 


16.0 


0.00541489 


2.6 


3 


4.320 


14.8 


0.0136554 


5.7 


5 


3.633 


13.5 


0.0336807 


11.6 


7 


2.957 


9.8 


0.249094 


35.2 


10 


1.394 


5.8 


0.849977 


19.9 


20 


0.400 



the constraints imposed on large-scale structure model by 
void statistics comes mainly from rare voids. 

5.2 Void mass density profiles 

In Fig.(|2l we show the mean density profile using expres- 
sion l(45|l for R = 10/i~^Mpc and So = —0.9, and for 
R = 8fe"^Mpc and So = -0.8667 

In Fig.(|3J we present the most probable profiles for the 
cases mentioned above including the 68 % confidence levels 
for both profiles, this levels define a quite narrow region up 
to over 13/i~^Mpc. 

Comparing these two figures it is apparent that, al- 
though both profiles are very similar within the voids, the 
mean profile is substantially steeper than the most proba- 
ble profile at the boundary of the voids. This is due to the 
asymmetry between the upper and lower one sigma profiles. 



Note that the profiles given here correspond to an av- 
erage over all empty spheres with quoted So and R, while 
those in the simulations correspond to the largest empty 
sphere with the same So and R. This implies that the latter 
profiles should be somewhat steeper than the former ones for 
r values slightly larger than R (within the sphere, and for 
r values substantially larger than R they should be equal). 
It is not difficult to account for this effect, but we shall not 
consider it here since it is not very relevant^ 

C omparing with simulations (fig. 3 in llCottlober et alJ 
l2003li l we find them to be in very good agreement. In par- 
ticular, for the R — 10ft~^Mpc void, the fiatness of the pro- 
file within the void with a gentle descent toward the center 
(5(0) = —0.93) is found in our results, as well as the steep 
rise at the boundary. This good agreement strongly suggest 
that the spherical collapse describes correctly the dynamics 
of individual rare voids even when their inner underdensity 
is quite low. That the spherical collapse model may give such 
good results in several cases, like the present one, where the 
degree of spherical symmetry is not that high and the tidal 
field due to the outside matter is not negligible it is an in- 
triguing fact that may be explained by the cancellation of 
the aspherical effects due to the local matter distribution 
and the tidal field generated by distant matter (Betancort- 
Rijo 2004). 

It may be checked that for these profiles, both for the 
most probable one and for the confidence limits, the La- 
grangian radius: 

g(r) = r(l + <5(r))'''' (62) 

is a monotonically increasing function of r. Thus, shell- 
crossing does not take place and expression 1391 applies, so 
that our procedure is consistent. One could think that this 
implies that, at least for 68 % of the profiles, shell-crossing 
does not take place. This is very nearly true, but it must 
be observed that, in principle, profiles within the limits may 
have wiggles, so that shell-crossing could be likely to have 
taken place; although even in this case it will not be very 
relevant, in the sense that eg. 1391 still very nearly applies, 
for values of r where the confidence region is narrow. 

ft must also be noted that it is not strictly true that 68 
% of the profiles are contained within the 68% confidence 
region. This is merely the region generated by the confidence 
intervals for 5 at a fixed value of r as r changes. That is, at a 
given value of r, 68% of the profiles must be within the region 
although the fraction of profiles that lay within this region 
for all the r values considered may be somewhat smaller. 
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Figure 4. The most probable enclosed density profiles for voids. In tiie left we sliow tlie profile for R = 8/i~^Mpc and (5o = —0.86667 
and in tlie riglit the plot for R = 10/i~^Mpc and Sq = —0.9. The solid line, for both, corresponds to the most probable profiles and the 
dashed line to the 68 % confidence levels. The dotted lines are the profiles found for voids in the numerical simulations. 




5 10 15 20 25 30 



Figure 3. The mean enclosed density profiles for voids. The thin 
line corresponds R = 8/i~^Mpc and Sq = —0.8667 while the thick 
line to R= lOh-^Mpc and Sq = -0.9 

5.3 Halo number density profiles 

We have calculated the most probable local halo number 
density profile, 5'j^f{r), within voids with R = 8/i~^Mpc 
and So = —0.8667, for masses above lO^/i~^M0 and above 
2 X lO^''/i~^M0, which correspond approximately to ha- 
los with circular velocity 20km /s ^ Vdr ^ 55km/ s and 
55km/ s ^ Vcir ^ 120km/ s respectively. To obtain 5jv(7') we 
have used Eq.(56) with S{r) given by the most probable mass 
density profile (Eq.(53)) corresponding to the mentioned val- 
ues of R and So- S'ff{r) has been obtained from (5jv(r) by 
means of Eq.(61). In Figure we show the halo number 
density profiles in and around the void and compared them 



with the local mass density most probable profile (obtained 
form 5{r) using expression 15811 ). It is apparent that \5n\ is 
larger than \S\ and the more so the larger the mean mass of 
the halos in the sample. In Figure JHJ we show the profiles 
presented in Fig.JSJ, but averagi ng it over five bins of e qual 
volume. Numerical results from iGottlober et alj (l2003l) are 
also shown for comparison. Although this last results corre- 
sponds to a superposition of five different values of R and 
So the agreement is quite good. Note that using our CMF 
is essential to explain the results found in the simulations: 
the ratio between the local halo number density (for masses 
above 2 x 10^'^ h~^MQ) at the center and at the boundary 
is about 2.6, while for the mean mass density this ratio is 
only 1.54 (for R = 8/i"^Mpc and So = -0.867). The ex- 
tra factor 1.69 is due to the different statistical clustering 
of the protohalos at the center and at the boundary, that 
is, to the dependence on position of the local number den- 
sity of protohalos before mass motion (i.e. on Lagrangian 
coordinates) . 

In figure Q we show the halo mass function for two dif- 
ferent voids. Note the excellent agreement with simulations 
(fig. (5) in Gottlober et al. 2003). 

Summarizing, the distribution of matter and halos of 
given masses within and around voids in simulations may 
be both reproduced by the combined use of the spherical 
expansion model and our CMF expression. 



6 DISCUSSION 

So far in this work we have been dealing only with voids 
defined by dark matter halos. However, the number density 
of voids defined by galaxies may be obtained in the same 
way as those defined by DM halos. 

To obtain the number density of the latter we implicitly 
had to determine the relative biasing of halos above certain 
mass with respect to the matter. This biasing was responsi- 
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Log M [Mg/h] Log M [Mg/h] 

Figure 7. Mass functions of halos averaged within two different voids. Left Panel: The thick full line corresponds to voids with radius 
R = 10/i-lMpc and mean density, &o = —0.9 obtained with our formalism. Also in this panel we show the mass function for 2 voids of 
the same parameters (R,5q) obtained from numerical simulations iGottlobcr et al. 2003)(thick dashed line). Right Panel: The thin line 
denote the mass function for voids with radius R = 8/i~^Mpc and 5q = —0.86667. Again, we show results for 3 voids with the same 
parameters from the numerical simulations. We can see that the agreement of our results with the numerical simulations is very good. 
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Figure 5. In this plot we show the fractional local halo number 
density within and around a void with R = 8/i~^Mpc and So = 
—0.867. The dashed line corresponds to halos with mass above 
10^/i~^Mq and the dotted line to those with mass above 2 X 
IO^^/i^^Mq. As a comparison we also show the local mass density 
profile (full line). 

ble for the fact that instead of using in expression (5jv = 5 
which correspond to objects distributed like mass, we had 
to use: 

1 + Sn = il+S){l+Sns) (63) 

where l + Sns, which is due to the initial statistical clustering 
of the protohalos, accounts for, or rather, is the origin of the 
biasing of halos with respect to mass. 

1 + S„3 was obtained by studying the dependence of the 
Lagrangian fractional fluctuation of the number of proto- 



.0 2 .4 .6 OB 1.0 

R/R 

Figure 6. The same as Fig. 4, where < 1 + <5^ > is the mean 
halo number density within spherical shells with 1/5 of the vol- 
ume of the void. The thin line correspond to halos with mass 
above 10®/!,~^Mq and the thick line to halos with mass above 
2 X lQi^^h~^yiQ. R/ R^aid is the distance to the center of the void 
in units of void radius. The open and half-filled squares corre- 
spond to halos with mass above 10^ and 2 X 10^^h~^yiQ respec- 
tively obtained for 5 voids by numerical simulations (Gottlober 
et al. 2003). The error bars denote the sampling error. 



halos within a sphere on the linear density fluctuation, 5i, 
within it. 

For voids defined by galaxies of certain type above a 
given luminosity, L, we must use in expression 16311 1 + 5ls, 
which describes the initial statistical clustering of the proto- 
galaxies, instead of 1 -I- 5ns ■ 1 + 5ls is obtained by means of 
the unconditional, nu{L), and conditional, nc{L, q, Q, 5i), lu- 
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minosity function in the same way as 1 + 5ns was derived by 
means of the conditional and unconditional mass functions. 

To obtain nu{L) and nc{L), we shall first assume that 
there exist a universal (independent of t he environment] 
Conditional Luminosity Function, 0(L/m) iMo et alJl2004l 
CLF). Then we have: 

poo 

nu{L) = I (j){L/m) nu{m) dm (64) 

nc{L,Q,q,5i) ^ / (j){L/m) nc{m,,Q,q,5i) dm (65) 
Jo 

Integrating over luminosity from L to infinity we obtain 
nu{> L) and nc(> L,Q,q,Si). Dividing the latter by the 
former we obtain 1 + SLs{q,Q,Si). Averaging over q within 
the void we finally have: 

f-i , X fr^xw 3 f'^ nc{> L,Q,q,5i) 2 , .c^n 

{l + 5LsiQ,5i)) = -— q dq 66 

Q Jo nu{> L) 

now, since ric depends on q and Q almost entirely through 
q/Q, in practice, 5ls depends only on Si. 

So, the number density of voids with radius larger than 
r defined by galaxies with luminosity larger than L is given 
by expression Q with Po{r) given by eq. @ and 1 + 5n 
given by: 

1 + SNiSi) = (1 + D{Si)){l + Sls{Si)) (67) 

More generally, one might consider the plausible possi- 
bility that void galaxies are a systematically different pop- 
ulation (Szomoru et al. 1996a; El-Ad & Piran 2000; Peebles 
2001; Rojas et al. 2004) or, in mathematical terms, that the 
Conditional Luminosity Function depends on the environ- 
ment. Evidence to this dependence on theoretical ground 
have been recently pointed out by Gao et al.(2005) who 
reported a dependence of halo clustering on environment 
through the environmental dependence on halo formation 
history. 

Assuming this dependence enters only through the envi- 
ronmental density we should then use (j){L/m, 82), where ^2 
stands for the present linear value of the fractional density 
fluctuation within a sphere centered at the galaxy and with 
radius 3 — 4/i~'^Mpc (this radius should be large enough to 
define a local environment but substantially smaller than the 
scale length on which these environment change, i.e. large 
voids) . 

The only change with respect to the previous case is 
that now to obtain nu{L) and nc{L) expression (25) must 
be multiplied by the probability distribution for S2 at q, 
P{S2/q,Q,Si), and integrated over 82- 

Up to here we have been using a non-uniform Poissonian 
clustering model: the probability per time unit for a galaxy 
to form at a given halo may be a function of time and some 
underlying field, but does not depend on whether or not 
some galaxies has previously formed in its neighborhood. 
Dark matter halos formed from Gaussian initial conditions 
may be shown to obey a Poissonian model. This is a conse- 
quence of the validity for th em of the peak-background split- 
ting approximation diBardeen et al.. .l986. . BBKS). However 
for galaxies this does not need to be true if the conditional 
luminosity function depends on the presence of neighbor- 
ing galaxies: A general model containing these possibilities 



is given in iBetancort-Riiol (I2OOOI) . To work within this gen- 
eral model we only need to change exp(— nV) in eq.Q by 
(1 + wn)~^^^ where w is an additional parameter to be de- 
termined from observations (for Poissonian model w = 0). 

By means of these expressions we may be able to use 
void statistics to impose constraints on the possible de- 
pendence of <j){L/m) on environment, and determine the w 
value. 



7 CONCLUSIONS 

The extension of the unconditional mass function that we 
have developed increases very substantially the reach of the 
original useful tool. Not only is the extension formally ex- 
act under the generally satisfied conditions that we have 
discussed, but it provides the local number density of col- 
lapsed objects at any distance from the point at which the 
constraint is evaluated. This is an essential point for the 
present work as well as for several other applications, since 
the fractional number density varies strongly from the cen- 
ter to the boundary of the void (the first value is roughly 
the cubic root of the second one for any mass). Previous pro- 
cedures for obtaining the conditional mass function assume 
that the mass function and the condition are evaluated at 
the same point. In principle, one could follow those proce- 
dures to obtain the mass function at any distance from the 
center of the void. But this is rather complex and can, at 
most, provide a good fitting formula. However, as we have 
shown, once the unconditional mass function has been ob- 
tained by deriving a fitting formula through the mentioned 
procedures (or any other) and calibrating it by means of 
numerical simulations, its extension to the conditional case 
can immediately be obtained, without having to repeat the 
derivation of the fitting formula and its calibration. 

This formalism has allowed us to obtain the mean frac- 
tional number density of collapsed objects of given mass 
within a void as a function of the fractional mass density 
within it. This has been used within the general formalism 
that we have described here to obtain the number density of 
voids with radius above a given value. We have compared our 
results with those found in numerical simulations, checking 
that for sufficiently rare voids our procedure gives very good 
results. Furthermore, using P{5/r) as given by the spheri- 
cal expansion/collapse approximation seems to be enough 
within present uncertainties. Note, however, that we have 
not checked separately to a sufficient extent the accuracy 
of the relationship between Po{r) and Po{r) (Eq.(l)) on the 
one hand and the accuracy of our computation of Po{r) on 
the other. We intend to do this by means of more detailed 
simulations that will allow us to eliminate some minor un- 
certainties thereby increasing the accuracy of our procedure. 

One relevant issue that we have not addressed here is 
the redshift distortions. Note that our formalism correspond 
to real space, while observations are made in redshift space. 
Recently, a procedure for correcting the observed statistics 
for redshift distortions has been developed (Patiri et al. 
2005). However, in a future work we intend to complement 
our formalism so that we could make predictions of voids 
statistics directly in redshift space. 

Our formalism provides a simple relation between prop- 
erties of the galaxy formation process and galaxy distribu- 
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tion statistics. Using this relation we may infer those proper- 
ties from the observed distribution. In this manner, in order 
to asses the consistency of a model of galaxy formation with 
void statistics it is not necessary to carry out complex nu- 
merical simulations with a large dynamical range but only to 
demand those models to show the properties (i.e. the condi- 
tional Luminosity function) inferred through our procedure 
from void statistics. Furthermore, the effect of the change of 
any parameter of the model may be estimated immediately. 
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APPENDIX A: ALTERNATIVE DERIVATION 

OF Sns 

In an alternative and more explicit procedure, instead of 
fitting directly the dependence of 5ns on m and 5i , we first 
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obtain the dependence of the fractional halo number density 
at q = 0, which we call S'^s, on m and Si. That is, before deal- 
ing with the mean fractional fluctuation within the sphere of 
radius Q we deal with its value at the center of this sphere. 
We find that 5'^^ may be approximated very accurately (for 
5i < 0) by: 

1 + C(5i,m) = ^'(771)6-*''™)*' (Al) 
where for m between IO^/i'^Mq and 2 x lO^^/i'^M©: 
1 

0.5 



A'(:m) 

b'im) = 0.0205 + 0.1155 



^3.51 X 10"/i~^Mq 
for larger masses we have: 
A'{m) ~ 1 

6'(m) = 0.1917 + 0.0198 



(A2) 
(A3) 



3.51 X 10"/i"^Mq 



(A4) 
(A5) 



To obtain the fractional density fiuctuation of the num- 
ber density of collapsed objects with masses above m, at a 
distance q from the center of the sphere under consideration, 
Snsim, 5i, q) we simply need to note that: 



o-i2(m, 0) 

<712(m.<;) 



(A6) 



where in the first parenthesis we have noted that in the 
procedure for extending the UMF expression the dependence 
on 5i enters essentially through 



ai2{m,q) 



Si 



(A7) 



which as we have seen (ea l35ll is independent of m. The 
second parenthesis accounts for the differences between the 
expressions that substitute d\na/dm in the extensions of 
UMF expression at g = and at a given q. This parenthesis 
is, for the values we shall use, very close to 1 and may be 
neglected. Having Snsim, Si, q) at any q we may immediately 
obtain Sns{m,Si) by computing its mean value within the 
sphere: 



APPENDIX B: DERIVATION OF P(5/i?, 5i) 

To obtain expression 1421 for the probability distribution for 
the mean value of S within r, we compute first the probabil- 
ity that S ^ A (where A is some given value) , which we rep- 
resent by Pc(A). If at a certain q value, go, Si{qo) (the linear 
value of the fractional density fiuctuation within Lagrangian 
radius q) were equal to DL{A) (the linear value correspond- 
ing through spherical collapse to an actual value A), then 
the sphere concentric with the void and with Lagrangian ra- 
dius go show a mean inner fractional fiuctuation equal to A 
corresponding to an expansion by a factor (1 + A)"^''^ and, 
if its present radius is r, we must have go = r{l + A)^''"'. If 
5; (go) were less than DL{A) then the actual S within this 
sphere (with Lagrangian radius go) will be less than A. So, 
the sphere with Lagrangian radius, go, would have expanded 
by a factor larger than (l+A)"^''"^ its present size being then 
larger than r. 

Thus, the sphere with present radius r would have ex- 
panded from an initial one with g' < go (neglecting shell 
crossing) and, assuming that the profile is monotonically 
increasing (which is obviously true for the most probable 
profile, but is also so for any realization with non-negligible 
probability), Si{q') shall be less than Si{q), and consequently, 
S{r) shall be less than A. We may then write: 

Pe(A) = P(5^A,/r,5i) = P(5, ^M(A),/go = r(l+A)^/^5l)(Bl) 



P{S/r,Sr) = —P,{A) 



(B2) 



where (Bl) denote the fact that the probability that S be 
smaller than A at r is, by the above arguments, equal to 
the probability that Si be smaller than DL{A) at go. But 
for fixed go and with the condition S{R) — So, Si follows 
distribution 13911 . so, expression (43) follows immediately. 

We have assumed above that shell crossing does not 
occur, however this may be shown to be the case in the 
relevant applications. Note that the relevant shell crossing 
here is that on the scale of the void, that is the shell crossing 
experienced by the mass field filtered on the scale of the void. 
On much smaller scales there is obviously shell crossing, but 
this is irrelevant to our argument. 



1 + Sns{m, Si) 



A (m)e ^ ' ' 



(A8) 



here, we have used eg. 11811 and eq.(Al), setting u = 
q/Q. We find again that this expression may be accurately 
fitted by ea. lAH with A{m) ~ A'{m) ~ 1 and 6(m) approx- 
imately equal to 6'(m)/2 in all the range of masses consid- 
ered. Sns is the mean (within Q) fractional fluctuation of the 
number density in Lagrangian coordinates (i.e., before the 
halos move along with mass). This fiuctuation is due to the 
statistical clustering, in the initial field, of the protohalos 
under consideration (those with mass above m). 
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